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ABSTRACT 

Among the many intriguing aspects of optically discovered tidal disruption events 
(TDEs) is that their temperatures are lower than predicted and that the temperature 
does not evolve as rapidly with decreasing fallback rate as would be expected in standard 
disk theory. We show that this can be explained qualitatively using an idea proposed 
by Laor &: Davis in the context of normal active galactic nuclei: that larger accretion 
rates imply stronger winds and thus that the accretion rate through the inner disk 
only depends weakly on the inflow rate at the outer edge of the disk. We also show 
that reasonable quantitative agreement with data requires that, as has been suggested 
in recent papers, the characteristic radius of the tidal stream is approximately equal 
to the semimajor axis of the most bound orbit of the debris rather than twice the 
pericenter distance, which would be expected from circularization without rapid angular 
momentum redistribution. If this explanation is correct, it suggests that the evolution 
of TDEs may test both non-standard disk theory and the details of the interactions of 
the tidal stream. 

Subject headings: accretion, accretion disks — black hole physics — galaxies: nuclei 


INTRODUCTION 


Tidal disruption events (TDEs), in which a star is torn apart by a massive black hole, have 
emerged in the last few years as excellent examples of what can be learned in time domain as¬ 
tronomy. The pos sible scient ific returns include elucidation of the powering of some active galac- 
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2011: Gezari et al. 
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2014 


Holoien et al. 2014 


Vinko et al 


20151 ) surveys. Among the many puz¬ 


zles posed by optically discovered TDEs is that their temperatures are not as high as predicted, and 
that the temperatures do not appear to decrease as rapidly as e xpected in standard disk th eory. 
For example, if a black hole of mass M accretes at a rate M, then IShakura &: SunvaevI (jl973l l find 
that the surface radiative flux at a radius R is 


^ ^ SGMM 

^ = “ 8 ^ 5 “ 


1 - 


Rin 


1 / 2 ' 


( 1 ) 


Here a is the Stefan-Boltzmann constant, G is Newton’s constant, and Rin is the inner radius 
of the disk, which in this formulation is assumed to be torque-free. Maximizing the flux over 
R we find that the peak disk effective temperature Tmax at a fraction / of the Eddington rate 
A^Edd = 4:7rGM/{Kr]c) is given by 
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( 2 ) 


Here n is the opacity, rj = L/Mc^ is the accretion efficiency that gives a luminosity L for an accretion 
rate M, and c is the speed of light. If we set R-^ = diGM/c^, which is the Boyer-Lindquist radius 
of the innermost stable circular orbit for a nonrotating black hole, then 


Tmax = 4.8 X 10® K /^/^(k/ 0.4 cm^ g-i)-i/^(77/0.057)"^/^(M/10® 


(3) 


(for a similar expression see, e.g., IUlmeiill999l . who also gives a somewhat lower temperature and 
different scaling if Rin is set by the tidal radius). Here we have scaled to the pure hydrogen 
Thomson opacity and the efficiency rj = 0.057 for a nonrotating black hole. Thus one would expect 
maximum disk temperatures Tmax ~ few x 10® K, and one would also expect that the temperature 
would decrease significantly over the course of an outburst. 

Neither of these expectations is met in optically discovered TDEs. Th e typical best-fit Planck 


temperature to the spectra is few x 10^ K rather than few x 10® K ( van Velzen et al 


2011 


Gezari et al.ll2012l : lArcavi et al.ll20l4 IChornock et al.ll20l4 iHoloien et al.ll2014l : IVinko et al.ll2015l l . 
It has been proposed that when the fallback rate is super-Eddington there is an optically thick 
photosph ere at many times the tidal radius, and that this could lead to characteristic temperatures 
~ lO'^ I£J Loeb fc Ulmer 1997 ; Ulmer et al. 1998 ; Strubbe fc Quataert 2009 ; Coughlin fc Beeelman 


20141 : ICuillochon et al.l 12014 ): see also Figure 3 of Stone fc Metzgei ( 20141 ) for a comparison with 
data. However, there are now multiple events in which the best-fit temperature remains roughly 
constant even as the fallback rate decreases by a large factor. For example, PSl~10jh maintains 
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a nearly constant UV-optical spectral shape even a s the integrated flu x from this event drops by 


more than an order of magnitude (see Figure 3 of iGezari et alJl2012l l. Similarly, ASASSN-14ae 
maintains a best-fit Planck temperature of ~ 2 x 10^ K (uncorrected for e xtinction) as the inte¬ 


grated flux drops by one and a half orders of magnitude (jHoloien et alJl2014l l. Standard disruption 


ture of the star in Phinnev 

1989 and 

Lodato et al. 2009 

and the relation between flux in a band 

and bolometric flux in Istrubbe &; Ouataert 

2009 and 

Lodato k Rossi 20111 or even as raoidlv as 

t if the stellar core survives ( 

Guillochon & Ramirez-Ruiz 

20131 suggest that the fallback rate 


would have dropped below Eddington after such a large drop in flux for black hole masses greater 
than ~ 10® Mq. It is possible that the size of a repro cessing layer adjusts wi th the fallback rate 
to yield a nearly constant maximum temperature (e.g., iGuillochon et alJl2014l ). but it is not clear 
what physics would enforce this. Thus explanations that involve super-Eddington accretion rates 
and expanded, nearly spherical, photospheres might not suffice for the entire observed evolution 
for several of these sources. 

Here we propose that both the lo w temperature and th e slow evolution of temperature can be 
explained in a framework proposed bv iLaor &: DavisI (j2014l ) to explain why the maximum temper¬ 
ature seen in AGNs is also only few x 10^ K over a wide range of black hole masses and accretion 
rates. The idea is that line-driven winds, which iLaor &: DavisI (120141 1 calibrate using data from 
stars, remove mass from the disk. At larger accretion rates at the outer disk boundary, more mass 
is removed and thus the peak temperature of the disk is reduced substantially from its standard 
disk value. At smaller accretion rates, less mass is removed and the temperature reduction is not 
as great. Thus the maximum achievable temperature is signihcantly less than it is in the standard 
disk picture, and the dependence of the peak disk temperature on the accretion rate is also much 
weaker than in the usual model. 

An additional quantitative test of TDE scenarios is available through the optical and near ultra¬ 
violet spectra that have been obtained for some events. If a star on a nearly parabolic orbit is tidally 
disrupted at a pericenter distance Vp and the debris circularizes at the specific angular momentum 
of the initial orbit, t hen the Newtoriian ci r cularization radius is 2rp. In contr a st, as is shown in sev¬ 
eral recent papers (IBonnerot et al.l 1201,^ : IGuillochon &: Ramirez-R.uizl 120151 : iHavasaki et al.l 120151 : 
Shiokawa et al.ll2015h if shock cool i ng is inefficient then the characteristic disk radius is much 
larger. For example, IShiokawa et al.l (1201.51 ) specifically suggest that the radius is roughly equal to 
the semimajor axis Omin of the most bound orbit. Given that Omin ~ (Af/m)^/^rtide) w here m is th e 
mass of the star prior to disruption and rtide for a stellar radius iZ* (|Reeslll988ll . 

the two radii can differ by a factor of tens if M/m ~ 10® — 10®. As we show, this has significant 
implications for the optical and ultraviolet fluxes, although precise predictions are susceptible to 
uncertainties in the shape and surface density distribution of the disk as well as to the unknown 
extinction in the host galaxies of TDEs. 


In Section [2] we review the salient aspects of the iLaor fc DavisI (|2014l ) model and show the 
expected evolution of the maximum disk temperature as a function of accretion rate. In Section [3]we 
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perfo rm by-eye fits to data from the PSl-lOjh (jCezari et alJl2012l ) and PSl-llaf (jChornock et al 
201J) candidate TDEs, which show that in the context of this model a large disk radius is favored 
strongly. We present our conclusions in Section HI 


2. METHODS 


2.1. Physical Principles 


Our suggestion is based on the disk wind model that iLaor &: DavisI (j2014l i propose to explain 
the few X 10^ K maximum disk temperature seen from AGN. Their wind loss prescriptions originate 
from observations of O stars, from which they derive two possible relations for the wind mass loss 
rate per area S. The first is: 


S = 2.6 X 10 g cm ^ s ^ 


(F/Fe) 


1.9 


(4) 


where F is the energy flux from a part of the disk and Fq = 6.3 x 10^^ erg s ^ cm ^ is the solar 
radiation flux. We refer to this as S(E). The second relation is 


S = 2.48 X g cm“2 s“^ 


{F/F^y-^^ {g/g^) 


- 1.11 


(5) 


where g is the surface gravity at the photosphere of part of the disk and Qca = 2.74 x 10^ c m s 


-2 


is the solar surface gravity. We refer to this prescription as 'E,{F,g). As iLaor &: DavisI (j2014l ) 


show, the assumptions of hydrostatic equilibrium and opacity dominated by electron scatter¬ 
ing imply that g/gQ = 2.5 x 10“^ F/Fq, which means that i^{F,g) reduces to S = 3.05 x 
10“® g cm“^ s“^ (F/Eq)^'^^. The second relation is a better fit to O star data, but it is not 
clear which, if either, is applicable to accretion disks around massive black holes. It is important 
to note that the power-law relation between S and F cannot, of course, extend to unlimited fluxes 
and effective temperature s. The true rela t ion de pends on the opac ity as a functi o n of t emperature. 
For example. Figure 1 of Cantiello et al. ('2OO9I) and Figure 38 of Paxton et al. ( 2013 1 display as¬ 
pects of the OPAL opacities ([iglesias &: Rogerslll996l l and show that above ~ 2 X 10^ K the 
opacity actually drops due to weakening of lines caused by high ionization. Thus at sufficiently 
high ionization fractions (which could be due to locally high temperatures or, as is probably the 
case in stellar-mass X-ray binaries, reprocessed X-ray emission from the inner disk), the wind loss 
rates will be suppressed. 

Because the winds fro m O stars are driven largely by the opacity due to the overlap of atomic 
lines (see the discussion in iLaor &: DavisI l2014l L there should be a significant dependence of the 
wind strength on the metallicity of the gas, with higher metallicity imply ing stronger w i nds i f 
other parameters are kept constant. This is not explicitly included in the iLaor &: DavisI (j2014l l 
model, but the theory could be adjusted if there is enough information about the metallicity of 
a particular disrupted star. We note that because metal line overlap is needed, the presence or 
absence of hydrogen is not expected to make a significant difference other than by allowing the 
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heavier elements to be a larger fraction of the overall ma ss. Thus the winds from the disruption 
of helium-rich stars, as has been proposed for PSl-lOjh bv iGezari et al.l (|2012lb would be expected 
to be stronger than but qualitatively similar to the winds from the disruption of stars with solar 
composition. 

Another point to consider is whether the wind will be optically thick and, if so, how it will 
interact with the outgoing radiation. To assess this we first need to determine whether atomic 
opacities or Thomson scattering dominate. A useful parameter in t his context is th e ionization 
parameter ^ = dvri^/ne, where Ug is the electron number density (see IDavidsonl 119771 for an early 
application to quasar s pectra). This is a measure of the number of ionizing photons per electron. 
Matt et al.l ( 19931 . Il996l l found that when ^ ^ 100 erg cm s“^ the gas is only weakly ionized, whereas 
when ^ ^ 5000 erg cm s“^ the gas is strongly ionized, which means that Thomson scattering (in 
this case, rather than free-free absorption) will be the most important source of opacity. 

If the radial speed Vr{R) ~ {GM/, then the number density of electrons in a spherical 
wind composed of pure hydrogen is roughly 


He. = 


M/rup 


M/nip 


4.'KR?Vr{R) 47ri?2(GM/i?)V2 ' 

The ionization parameter is ^ = AirF/ne = At:{L/A'KR?‘) ne-, so if we use L = r]McP we find 
^ = rimpC^A7r{GM/^ 2 x 10® erg cm GM/c^)/R]^^^ . 


( 6 ) 


(7) 


Thus the wind is highl y ionized and Thom s on sc attering is the dominant source of opacity. In 
contrast, if we use the IShakura Sz SunvaevI (jl973l l solution for the inner region at R ^ i?in, we 
find ^ ~ 200 erg cm s“^ at the photosphere of the disk for M = 3 x 10® Mq, R = lO^GM/c^ , and 
M = Mejcj- This means that line opacities can dominate radiation transport in the disk but are 
unlikely to be important in the wind. 


The optical depth through the wind is then 


r = nardr = 

^ (0.06/r?)(M/MEdd)[(10® GM/c'^)/R]^/^ 


( 8 ) 


independent of the black hole mass. The numerical coefficient is roughly unity for r] = 0.057, which 
means that for our standard parameters the optical depth could be a few for mass fallback rates of 
~ lOx Eddington. However, the very high ionization parameter implies that absorption opacities 
will not be important in the wind. Thus the spectrum might be mildly Comptonized but seems 
unlikely to be affected otherwise. However, as the referee pointed out, this scattering could smooth 
out variations in the light curve that are shorter than the diffusion time through the wind, when the 
wind is optically thick. This could be hours to days in the peak fallback rate part of the disruption. 
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2.2. Numerical Implementation 


For simplicity we assume that the disk has Keplerian rotation, i.e., = GM/R^ at a disk 

radius R. This requires that the disk is geometrically thin, because otherwise pressure support 
will change the rotation law. The thin-disk assumption will be poor close to the hole if the local 
accretion rate is comparable t o or la rger than Eddington, but at hundreds of gravitational radii the 
standard IShakura fc SnnvaevI (j 197, 'll ) solntion snggests that even disks with accretion rates several 
times Eddington will be geometrically thin. 


With this assumption, manipnlation of equations (26)-(28) of iLaor &: DavisI (j2014l i yields the 
differential equations to be integrated: 




aw, 


dR 


Rtp 


M 

AttR^ 


m 


1/2 


-m 


(9) 


R<t> 


where is the vertically integrated stress and as above S is a fnnction of F, where 


F = -W, 


R<t> 


eGM^ 

\1f) 


2 R 


( 10 ) 


Here e, a parameter in the iLaor &: DavisI (j2014l l model, is the ratio of the wind speed at infinity 
to the escape speed from the location in the disk where the wind originated. Note that although 
we treat e as independent of the effective temperatnre, it need not be: for example, as pointed out 
by Nick Stone (private communication), there is some evidence fro m massive stars that the re is a 


jnmp in wind speed when the temperature is ~ 1 — 2 x 10^ K (e.g., IPauldrach fc Pulslll99fll ). 


Solution of this system of equations requires a choice for the inner boundary condition on 
Eor simplicity we choose the zero-stress zero-spin bonndary condition 


= 0 at i?in = GGM/ 


( 11 ) 


and verify that when M is set to zero everywhere the solntion agrees with the analytical expression 
Wr^ = (MHfc/27r) [l — (i?in/7?)^'^^]. Note that zero inner str ess is neither phy sically reqnired 
nor consistent with all magnetohydrodynamic simulations (e.g., iNoble et al.l 12011)1 ): it is simply a 
convenient reference case, and the resnlts we display depend only weakly on the inner bonndary 
condition. Eor numerical reasons it is easiest to integrate outward from the inner boundary rather 
than inward from the outer boundary. Thus for a given desired mass accretion rate Mout at the 
outer disk boundary Rout we make a guess for M\^ at the inner boundary R\a (which we set to 
QGM/c^ for the resnlts shown here), integrate outward to i?out and measnre Moutj and converge 
on the solution for Min by bisection. 

In Eigure m we show the effect of the wind for various prescriptions (E(E) or T,(F,g), and 
e = 0 or e = 1). We assume that the outer accretion rate equals the Eddington rate, and track the 
luminosity gradient times the radius, which shows the radius at which the luminosity peaks. We see 
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that at large distances from the hole the solutions converge, as expected. At smaller distances, the 
luminosity is less at a given radius than it would be without a wind, by an amount that depends 
somewhat on the prescription. We also see that the luminosity peaks at a larger radius than it 
would without a wind. One consequence of this is that the maximum temperature of the disk is 
lowered significantly by a wind. 


This maximum temperature is plotted in Figure [2j We show the disk maximum temperature 
as a function of outer accretion rate, assuming S = T,(F) and e = 0; other choices for S and 
e give qualitatively similar results. The particular dependence we find is close to the g s 

6.5 X K [{M/M^dd)/{M/10^ dependence found analytically by Laor &: Davis ( 2014 1. 

We see that wind losses weaken the dependence of maximum disk temperature as a function of 
outer accre tion rate, as is consistent with several r e cent observations o f optically detected TDE 
candidates (jGezari et al.ll2ni2l : IChornock et al.l 120141 : iHoloien et al.ll2^14^ . 


RESULTS 


In this section we display illustrative hts of our model to data from two optically detected 
TDEs. In Eigure[3l we sho w an exampl e comp arison of the E(E), e = 0 model to data from the 
PSl-lOjh TDE taken from iGezari et al.l (j2012l l. Here we assume a nonrotating black hole with a 
mass of M = 3 x 10® Mq (other masses give comparably good fits), and we have incorporated 
the redshift z = 0.1696 and the implied luminosity distance d h,rr , = 820 Mpc to this source. The 
dotted line shows the disk wind model spectrum for a disk with an outer radius comparable to 
the circularization radius for a stream with a pericenter equal to the tidal radius. For this line we 
assume an accretion rate at the outer edge that is 15 times the Eddington accretion rate for a black 
hole of this mass that accretes pure hydrogen. Clearly the optical and ultraviolet fluxes are far too 
small, and the colors are far too blue, compared with the data. In contrast, the solid lines show 
disk wind models with an outer radius Rout = 1400 GM/c^, which is comparable to the semimajor 
axis of the most bound debris. For the data 19.1 days before maximum, we assume an outer disk 
accretion rate that is 15 times the Eddington accretion rate for pure hydrogen, and for the data 
242.4 days after maximum, we assume an outer accretion rate that is 0.25 times the Eddington rate 
for pure hydrogen. These models are obviously far superior to the small-disk models, although the 
agreement with the data is not perfect. 

Similarly, we sho w in Figure IH a by-eye fit of the S(E), e = 0 model to data from the PSl- 


llaf TDE taken from IChornock et al.l (j2014l l. We again assume a nonrotating black hole with 


M = 3 X 10® Mq (and again, other masses would work as well), and we have incorporated the 
redshift z = 0.4046 and the implied luminosity distance dium = 2200 Mpc to this source. Eor the 
data at 8.2 days after maximum, we assume an outer disk accretion rate at the outer edge that is 14 
times the Eddington accretion rate for pure hydrogen. Eor the data from 66.4 days after maximum, 
we assume an outer disk accretion rate that is 1.8 times the Eddington accretion rate. In both cases, 
our outer disk radius is i?out = 3000GM/c^. The fits to PSl-lOjh and PSl-llaf suggest that the 






















heat from stream interactions cools inefficiently ( 

Bonnerot et al. 

1—1 

o 

CN 

Guillochon &: Ramirez-Ruiz 

2015; 

Havasaki et al. 

2015: 

Shiokawa et al. 

20151. 


The time scale of the evolution of the fallback rate implied by our fit to the PSl-l Ojh data is 
approximately consistent with models. For example, Figure 21 of IShiokawa et al.l (j2f)lf)l ) shows the 
ballistic fallback rate scaled from their simulations (this corresponds to the outer disk accretion 
rate in our model). The unit of time is the orbital time at the outer disk radius, which for our 
example fit is 2ti = 56 days for i?out = 1400 GM/(P and M = 3 x 10® Mq. The 
implied rise time of ~ 50 days is comparable to wh at iGezari et al.l (1201211 found for PSl-lOjh. The 


scale M^/fGM/ c^) for the accretio n rate used bv IShiokawa et al.l (120151 1 is 3.13 x 10^^ erg s ^ if 


M* = 0.23 Mq (jGezari et al.ll2012l l. so at 19 days before the peak, or ~ 0.3 time units before the 
peak, the expected accretion rate would be 2 x 10“® M*/(GM/c^) 6.26 x 10^® g s“^, which is 

roughly 10 times the hydrogen Eddington accretion rate onto a 3 x 10® Mq nonrotating black hole. 
At 242 days after the peak, or approximately 4.3 time units after the peak, the expected fallback 
rate would be roughly 10 times less than it was 19 days before the peak rather than the 60 times 
less that we find in our model. 

In Figure [5] we show an example of the outer and inner accretion rate inferred using a particular 
model of the PSl-lOjh event. Here we assume S = S(T), e = 0, M = 3 x 10® Mq, and Ront = 
1400 GM/(? for all data, and fit only t o the best-estimate g-band magnitude without extinction 


corrections given in iGezari et al.l (120121 ) . There are many reasons wh y the disk structu re could 


differ from our assumptions (for example, as was recently pointed out bv iPiran et al.ll2015l . the disk 
could remain elliptical for many orbits), but the general agreement with a Mout {t — to) law 
is encouraging. In this model the integrated mass that fell on the outer edge of the disk over the 
~ 1 yr of observation was Mout ~ 1.5 Mq, but the integrated mass that fell into the black hole was 
only ~ 2 X 10“^ Mq. Similar masses are inferred for S = or e = 1, and lower masses (Mout 

equal to a few tenths of a solar mass) are inferred for larger outer radii. The implied fallback mass 
is plausible for a stellar disruption, but the overwhelming majority of the matter is driven away 
in a wind. This wind will carry a considerable amount of kinetic energy Ekin ~ Mqu^ ~ 10®® erg; 
whether, and how, this energy becomes manifest will likely depend on the strength of interactions 
the wind has with any surrounding gas in the galactic center. 

The time scale implied by o ur fit to the PSl-llaf event does not fit as well as it does to the 
PSl-lOjh data. From the data of iGhornock et al.l (j2014l l the rise time of the event is shorter by tens 
of percent than it is for PSl-lOjh, but our model requires a larger disk and therefore a larger char¬ 
acteristic time. Perhaps in this source the self-intersection of the tidal stream proceeds in a slightly 
different way. For example, it could be that the shocks due to stream-self-inter section in PSl- 


llaf a re particularly strong; this drives the initial inflow, as was demonstrated bv IShiokawa et al 


(|2015l ). Thus for PSl-lOjh our model provides a reasonable fit to the spectra and to the time scales, 
whereas for PSl-llaf our model provides a reasonable fit to the spectra but is off by a moderate 
factor from the observed time scales. 
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Many caveats apply to such fits. For example, there is an unknown level of extinction intrinsic 
to the host galaxy, there is an unknown amount of contamination from other sources, and the 
smoothness of the ultraviolet and optical data means that many models give equally good formal 


fits (t his includes disk models without mass loss to winds and, as is shown in Figure 3 of iGezari et al 


2012I . even simple Planck spectra). There are also likely to be other sources of variation in mass 


accretion rate with radius, especially if the inflow time is long compared with the mass accumulation 
time at the outer edge of the disk. In addition, for this fit we assume for simplicity that each disk 
annulus emits as a bla ckbody, whereas in reality there will be spectral corrections introduced by 
radiative transfer (e.g., iDavis et al.ll2005l f. What we highlight here is what we believe to be a robust 
result; in order to explain the ultraviolet and optical data it is necessary that the characteristic 
disk size be considerably greater t han than the that comes from straightforward angular 

momentum conservation (see also IStone &: Metzgeu [201J) . This is true even relatively early in the 
disruption, which suggests that the disk becomes large quickly. 


CONCLUSIONS 


We have shown that two aspects of optically discovered TDEs can be explaine d qualitatively 
by a m odel involving accretion disk winds that was originally proposed for AGN bv iLaor fc Davis 
(j2ni4l ). These aspects are the low maximum temperature from TDE disks and the weak evolution 
of that temperature as a function of the fallback rate. In both cases the idea is that wind losses 
throughout the disk act to regulate the accretion rate through the inner portion of the disk. There¬ 
fore, increasing the accretion rate at the outer edge of the disk increases the inner disk temperature 
by a far smaller factor than would be expected in the standard no-wind disk model. This implies 
that massive black holes are messy eaters of stars, in the sense that the accretion rate is much 
lower than the supply rate from the disruption; for example, as we showed in Section [3l it could 
be that well under 1% of the infalling matter is actually accreted by the black hole. Thus con¬ 
sumption of stars may be a highly inefficient way of growing massive black holes. In addition, we 
find suggestive evidence in two TDEs that the characteristic disk size rapidly exceeds 2rp, which 
would be the expected circularization radius, and approaches the size Omin of the most bound 
orbit of the debris. This is consistent with recent simulations in which cooling of the shocked 


debris is inefficient (Bonnerot et al. 2015 : Guillochon fc Ramirez-Ruizi 


2015 


Havasaki et al 


2015 


Shiokawa et al 


Shiokawa et al.l 120151 ). and in particular consistent with the specific suggestion of 
(j2015ll that i?out ~ amin- One might, in this model, expect that the wind would have an effect on 


the profiles of broad emission lines (e.g., iFlohic et al.l 120121 show that such lines become narrower 
with increasing wind optical depth), but a detailed investigation of the implications for lines is 
beyond the scope of this paper. 


A natural question is: under what circumstances can the maximum disk temperature be high 
enough that it produces detectable X-rays, as is required for X-ray discovered TDEs? There are 
several possibilities, none of which is likely to suffice in every case. First, in some circumstances, it 
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could be that the matter is not sufficiently thermalized to emit as a blackbody, but instead radiates 
less efficiently and thus reaches a higher t emperature (as could happen in coronae or jets, or in 
the thermal bremsstrahlung suggestion of ISembav WestIIl993l l. Second, given that line-driven 
winds require heavy elements, if the metallicity is low then less mass is driven away and greater 
temperatures can be achieved. Third, as in stellar-mass black hole binaries, there could be a thresh¬ 
old temperature beyond which the ionization fraction is great enough that the relevant opacity is 
Thomson rather than the opacity given by overlapping metal lines; this would become progressively 
more important at lower black hole masses, so X-ray detection could sometimes indicate lower mass 
primaries. If this is the explanation then there could be circumstances in which, as the fallback 
rate drops and the temperature decreases, line opacities become more important and there is a 
rapid drop in X-ray flux as the disk temperature decreases. Finally, on occasion, variability of the 
accretion onto black holes could be mistaken for a TDE. 


In summary, disk winds offer a natural explanation for the low temperatures and weak de¬ 
pendence of temperature on mass accretion rate that is seen in several optically-discovered TDEs. 
If this explanation is correct, it suggests that observations of TDEs can probe non-standard disk 
models as well as the self-interaction of streams after the initial disruption. 
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Fig. 1.— Gradient of bolometric luminosity times radius as a function of radius for models with 
M = 3 X 10® Mq and an outer accretion rate equal to the Eddington rate. The solid line shows the 
standard no-wind solution, and the other line types show solutions with winds assuming S = S(F) 
or S = and e = 0 or 1. As expected, at large radii the solutions asymptote to each 

other, but at smal ler radii the mass los s reduces the luminosity considerably in the wind models. 
See Figures 4-6 of iLaor &: DavisI (|2014l i for similar comparisons of the accretion rate and effective 


temperature as a function of radius, and of the total spectrum from the disk. 
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Fig. 2.— Maximum disk effective temperatur e vs. infall accretion rate in Eddington units at the 
outer edge of the disk as predicted using the iLaor &: DavisI (j201J) model (solid line), with e = 0 
and S = S(E) (see Section [2] for details; other model choices give similar results). For comparison 
we show the maximum temperature in the standard no-wind disk solution (dotted line). Over a 
factor of 10^ in infall rate the maximum temperature in the wind model changes by only a factor 
of 1.9, compared with the factor of 10 predicted in the standard disk model. In addition, the 
maximum temperature at an Eddington accretion rate for this M = 3 x 10® Mq black hole is only 
7.4 X lO'^ K, rather than Tmax = 3.2_ X 10® K as predicted for standard disks. In both respects, 
the disk wind model of iLaor &: DavisI (120141 ) is a better fit than the no-wind model to data from 
optically discovered TDEs. 
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Fig. 3.— Illustrat ive fit of the disk w ind model to the magnitude data for the candidate TDE 
PSl-lOjh given in iGezari et al.l (|2ni2l ): the solid symbols with error bars are for 19.1 days before 
maximum luminosity, and the open symbols with error bars are for 242.4 days after maximum 
luminosity. See text for details of the model. The solid lines are for an outer disk radius that is 
comparable to the semimajor axis of the most bound debris, whereas the dotted line shows the 
output of the model assuming an outer disk radius that is twice the tidal radius. Ot her models fit 


equa lly well (e.g., a standard disk model with no wind, or even a blackbody as shown bv iGezari et al 


2012l i. but this hgure demonstrates that both the fluxes and the colors require that the emitting 


area must be closer to the area charac teristic of the most bo und debris than to the area at the 
nominal circularization radius (see also lStone Metzgeiil2014| j. 














21 


NUV 


g 


r i z 



Fig. 4.— Same as Figure[3l but for the candidate TDE PSl-llaf ( Chornock et al. 2014 ). The data 
are again labeled by time relative to maximum luminosity. For both lines the outer disk radius 
is 3000GM/c^. The outer accretion rate for the data 8.2 days after maximum is 14 times the 
Eddington rate for a pure hydrogen disk around a M = 3 x 10® Mq black hole, and for the data 
66.4 days after maximum is 1.8 times the Eddington rate. Again, the disk wind model provides 
reasonable fits to the spectral data. 
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Fig. 5.— Mass accretion rate at the outer edge of the disk (Mout, labeled with the solid line) and 
into the black hole (Min, labeled with the dashed line), for a particular model of the PSl-lOjh 
event in which S = S(F), e = 0, M = 3 x 10® Mq, and Rout = 1400 GM/c^ for all hts. We 
ob tained M by finding the accretion rate that gave a g-band magnitude equal to the best estimate 
in iGezari et al.l (j2012l '): beyond ~ 240 days the large photometric uncertainties lead to choppiness 
in this best estimate. There were no data from ~ 50 to ~ 180 days after the peak. For comparison 
we also plot an offset curve that shows the relative accretion rate for a (t — l^-w. Caveats 

apply regarding the simplicity of the model, but the inferred fallback rate is comparable to what 
is expected. 











